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Abstract 

In a previous paper we formulated a set of necessary conditions for the spherically 
symmetric weakly charged dust to avoid Big Bang/Big Crunch, shell crossing and 
permanent central singularities. However, we did not discuss the properties of the 
energy density, some of which are surprising and seem not to have been known up 
to now. A singularity of infinite energy density does exist - it is a point singularity 
situated on the world line of the center of symmetry. The condition that no mass 
shell collapses to ii = if it had R > initially thus turns out to be still insuffi- 
cient for avoiding a singularity. Moreover, at the singularity the energy density e is 
direction-dependent: e — > — oo when we approach the singular point along at = const 
hypersurface and e — > -|-oo when we approach that point along the center of symme- 
try. The appearance of negative-energy-density regions turns out to be inevitable. 
We discuss various aspects of this property of our configuration. We also show that 
a permanently pulsating configuration, with the period of pulsation independent of 
mass, is possible only if there exists a permanent central singularity. 

1 The problem 

In our previous paper [Ij (Paper I) we formulated a set of necessary conditions for the 
spherically symmetric weakly charged dust to avoid Big Bang/Big Crunch, shell crossing 
and permanent central singularities. In fact, it had already been proven by Ori [21 [3] that 
weakly charged dust with no central singularity necessarily evolves toward a shell crossing. 
However, Ori assumed that the absolute value of the charge density (in geometric units) 
is strictly smaller than the mass density throughout the volume, including the center. His 
proof did not include the subcase when the values of these two densities at the center of 
symmetry are equal. We discussed precisely this subcase. We found that shell crossings 
are still unavoidable when the energy function E > 0, but the conditions did not lead to a 
contradiction when E < 0. We also provided an example of a solution of these conditions, 
which was in fact incorrect (see below). However, we did not discuss the properties of 
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the energy density. While we started to investigate more properties of the numerical 
example presented in Paper I, it turned out that this density was becoming negative in a 
vicinity of the center of symmetry, for a brief period around the bounce instant. Then, 
exact consideration showed that this is a general problem - a weakly charged spherically 
symmetric dust ball will always have such a negative-energy-density region if the charge 
and mass densities become equal in absolute value at the center. This fact is hidden in the 
equations in a rather tricky way. In the general case treated by Ori this problem can be 
avoided by an appropriate choice of the arbitrary functions. 

The explicit example we gave in Paper I had a problem - it did not obey the condition 
that the ratio R/ Ai^^^ must have a finite nonzero limit at the center of symmetry. For our 
example, this limit was zero, which means that there was a permanent central singularity 
in it (the limit of the mass density at the center was actually — oo). The permanent infinity 
is easily cured by changing x^^^ to x^^^ in the formulae for P and E. The whole subsequent 
reasoning then applies with only little quantitative changesjl] but a direction-dependent 
point singularity at the center necessarily appears at the instant of maximal compression. 

In the present paper we fill in the gaps left out by previous investigation. We show that 
the energy density of the charged dust must become negative within a finite time-interval 
containing the instant of maximal compression. We also show that the point on the world 
line of the center of symmetry that is the limit of the points of maximal compression is 
a direction-dependent singularity of infinite energy density. The energy density tends to 
minus infinity when we approach the singular point along the hypersurface of maximal 
compression, and to plus infinity when we approach the same point along the center of 
symmetry. Thus, the condition that the areal radius R{t, r) never goes down to zero if it 
was nonzero initially turns out to be still insufficient for avoiding singularities. These are 
the most striking and most important results of the present paper. 

In order to make the paper independently readable, we quote the basic results of Paper 
I in the next section, without proof. All the proofs (or references to the proofs) and 
details of the calculation can be found in Paper I. Section [2] also contains remarks on the 
interpretation of the energy density, which is nontrivial in the presence of charges, and 
was not properly discussed in the earlier papers. In section [3] we show that the quantity 
u = r — QQ,fyf /i? must become negative in a vicinity of the instant of maximal compression 
whenever the conditions Q,n'^ = G/c^ at the center and Q^n"^ < G/c^ elsewhere are obeyed. 
In section H] we discuss the various consequences of this fact and show that this necessarily 
implies a negative energy density for a certain period around the bounce instant. In section 
E] we show that the negative energy density regions cannot be eliminated by assuming 
that the charged dust ball contains an empty (Minkowski) vacuole around the center of 
symmetry because such a vacuole cannot exist. In section [7] we prove the existence of 
the point-singularity on the world line of the center of symmetry. In section [9] we show 
that, even with negative energy density allowed, a solution with the period of oscillation 
independent of mass does not exist. A permanently nonsingular configuration of nonstatic 
weakly charged dust is thus impossible, unless one allows a permanent central singularity. 
We also compare our conclusion with the properties of the uncharged case. Section [TO] 
summarizes the conclusions. 

"'^This error has already been pointed out in an erratum to Paper I, and a fully corrected text is available 
from the gr-qc archive, see Ref. [I]. 
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2 Basic formulae and results of Paper I 



For a spherically symmetric spacetime in comoving coordinates, the metric can be put in 
the form 

ds^ = e^(*'")dt2 - e^(*'")dr2 - R^{t, r) [M^ + sm\^)dip^] . (2.1) 

In the generic case R,r 7^ 0, assuming there are no magnetic charges, the Einstein-Maxwell 
equations yield the following result [H IH [5] . The only independent nonzero component of 
the electromagnetic field is 

F^' = g(r)e-(^+^)/V/^^ (2.2) 
Q,r = {A7r/c)p,e^/'R\ (2.3) 

where Q{r) is an arbitrary function - the electric charge within the r-surface, and pe is the 
electric charge density; 

^ei?V/2 = ^iV„, /^ = 87^G/c^ (2.4) 

where e is the energy density and N,r is an arbitrary function of integration. The N so 
defined corresponds, in the electrically neutral case, to the energy equivalent to the sum 
of rest masses within the r-surface. With charges present, this interpretation involves a 
subtle point, see below after eq. (12.101) . The ratio Q,^ /A^,r = Pe/(ce) is time-independent. 

T{r)-QQ,N/R' ^^-^^ 



where Q^n is just an abbreviation for (5,r IN 



A/2 



C,r=2^QQ,N; (2.6) 

e-^R/ ^ - 1 + ^ + QHQ,.^-G/c^) 1 

R i?2 3 ' V ; 

where A is the cosmological constant, and M(r) is an arbitrary function. By analogy with 
the uncharged Lemaitre - Tolman model, (F^ — l)/2 is often denoted as E{r). 

The function M(r) is the effective mass, and it need not be positive. It is connected 
with the previously defined arbitrary functions by 

^FiV„=(M + QQ,^F)„. (2.8) 

The quantity 

M = M + QQ,N F, (2.9) 

is the active gravitational mass. Thus, via (12. Sp . F determines by how much Ai increases 
when a unit of rest mass is added to the source, i.e. F is a measure of the gravitational 
mass defect/excess. 

The energy density of the dust is: 

,, = if^fr-«^V (2.10) 



C^R^R,r V R 
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Now note the subtle point. In the electrically neutral case, Q = 0, the quantity e is the 
density of rest-mass (in energy units), so that the 3-space integral fy y/—gedsx simply 
equals the sum of rest masses in the volume V divided by c^. However, with charges 
present, e contains a contribution from the charges. As will be shown later in this paper, 
QQ,N = \ (Q^) ,T- iN^r must be positive in a vicinity of the center of symmetry, and R > 
everywhere for geometrical reasons. Thus, the presence of charges always decreases 
the energy density. It will turn out later that, in the class of models considered here, e 
necessarily becomes negative for a brief period around the instant of maximal compression 
along each world-line except the central one. 

The functions R{t,r) and C{t,r) are implicitly defined by the set (12. 6p - (12.71) . which 
can in general be solved only numerically. 

If the configuration considered here is matched to the Reissner - Nordstrom (R-N) 
metric across a hypersurface r = ri,, then the following must hold [U El [1]: 

^Qin), m = {M + QQ,^T)^^^^, (2.11) 



where e and m are the R-N charge and mass parameters, respectively. 

Unlike in the electrically neutral case, in charged dust the Big Bang/Big Crunch 
(BB/BC) singularity can be avoided, i.e. there exist solutions of (12. 6J) - (12.71) in which 



the function R{t, r) never goes down to zero if it was nonzero ini tiallyi We give the con- 
ditions for the existence of such solutions only in the case A = 0. One of those conditions 
is 

M^>2EQ^Q,,,^-G/c''), (2.12) 

which is fulfilled identically when E = 0. With (I2.12p fulfilled and E ^ 0, the right-hand 
side of (12. 7p has two roots, given by 

M 



The second condition for avoiding the BB/BC singularity depends on E: 

(a) When E < 0, the condition is 

Q,N^ < Gjc^ and M > 0. (2.14) 

With (I2.14P fulfilled and E < 0, R oscillates between the nonzero minimum at i? = R^ 
and the maximum at i? = -R_. With (12.140 and E = 0, R goes down from infinity to the 
finite minimal value -Rmin = ip/c^ — Q,n'^) /(2M) and then increases to infinity again. 

(b) When E > and M > 0, the BB/BC singularity does not occur if Q,^'^ < G/c^. 
When E > Q and M < 0, nonsingular solutions exist with no further conditions, provided 
R > i?+ initially. The bounce with M < is nonrelativistic, since it occurs also in Newton's 
theory, under the same conditions (M < means that the electric repulsion of the charges 
spread throughout the volume of the dust prevails over the gravitational attraction of the 
mass) . 



^'R{t,r) — permanently at the center of symmetry. 
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The inequality Q^n'^ < G/c^ translates into \pe\ < \/Ge/c, which means that, in geo- 
metric units, the absolute value of the charge density is smaller than the mass densityo 

There is another subtle point here. The conditions (12.141) guarantee that a particle 
that had i? > initially will not hit the set i? = in the future or in the past. But, as we 
will see in Sec. [TJ the configurations obeying (I2.14p contain a cleverly hidden singularity 
of a type hitherto unknown in dust solutions. On the world line of the center of symmetry, 
where R{t,r) = permanently, there is a point in which e +oo for a single instant. 
This instant is the limit at i? — > of the hypersurface Smin consisting of the instants in 
which the mass shells with R> attain their minimal sizes. However, if we approach the 
same spacetime location along the hypersurface Sram, then e — > —oo. 

The surface of the charged sphere obeys the equation of radial motion of a charged 
particle in the Reissner-Nordstrom spacetime. Thus, the surface of a collapsing sphere 
must continue to collapse until it crosses the inner R-N horizon r_, and can bounce at 
R < r_. Then, however, it cannot re-expand back into the same spacetime region from 
which it collapsed, as this would require motion backward in time. The surface would thus 
continue through the tunnel between the singularities and re-expand into another copy of 
the asymptotically flat region. 

At small charge density {Q,n'^ < G/c^ throughout the volume) a shell crossing is 
unavoidable, and it will block the passage through the R-N tunnel, as shown by Ori 
[2], 13]. A nonsingular bounce might be possible only if Q,n'^ > G/c^ everywhere or if 
Q,n'^ G/c^ at the center, while Q^n'^ < G/c^ elsewhere. The first case was dealt with 
by Ori (unpublished [6]). We consider the second case in Section 3, and the result is that a 
nonsingular bounce might happen only if the energy density becomes negative for a period 
around the bounce instant. 

Shell crossings are most conveniently discussed in the mass-curvature coordinates (A^, i?), 
first introduced by Ori [2]. Details of the transformation can be found in Ref. The 
{M.,R) coordinates allow the Einstein-Maxwell equations with A = to be explicitly in- 
tegrated, but it happens at a price. The spacetime points are identified by the values of 
M. and i?, which means a given point is defined by saying "it is the place where the shell 
containing the mass M. has the radius R\ The information on the time-dependence of 
R is lost, and can be regained only by reverting to the comoving coordinates - but the 
transformation equations are equivalent to (12.61) - (12.71) and cannot be exphcitly integrated. 

The solution of the Einstein - Maxwell equations is given below. The velocity field has 
only one contravariant component: 



(+ for expansion, — for collapse). We define the auxiliary quantities 

def „ ,„ . dcf . 2Ad GQ^ 1 A n2 /r^ 

u = T-QQ,n/R, ^ = ^"~^ + ^ + 3^^ ' (2-^^) 
and we get for the metric 

gMM = F^d., gMR = Fu/u'', (^^^ = l/(n^)2, (2.17) 



■^This includes also the case Q,n — 0, i.e. zero charge density, provided Q ^ 0, i.e. nonzero total charge. 
Such a configuration is neutral dust moving in the exterior electric field of a spherically symmetric source. 
Also in this case, the BB/BC singularity is avoided. 
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where the function F{M., R) is given by 

1 r r4 1-) 

(2.18) 



1 1 

\ r,_A/( +- 



c4 

-Q {Q^n'^ + QQ,NN ) 



With A = 0, the integral of this is elementary, but to give its explicit form several cases 
have to be considered separately (see Ref [2j for a list). Shell crossings occur at the zeros 
of the function F. Thus, to avoid shell crossings, the arbitrary functions must be chosen 
so that F ^ everywhere. 

The only nonvanishing components of the electromagnetic tensor in the {Ai,R) coor- 
dinates are 

pMR^_pRM^^^^ Fj^^ = -FnM = -^, (2.19) 

while the charge density and the energy-density are 

47rp, _ Q,M _ 2 

c ~ R^Fu^' TR^Fu^- ^ ' 

The (A1, i?)-coordinates cover only such a region where R,t has a constant sign. The 
function F changes sign where R,t does, and so does = e~'"/^i?,j. Thus, Fu^ preserves 
its sign when collapse turns to expansion and vice versa. 
In the (A^,i?) coordinates eq. (12. 8p reads 

M,N={GT/c''). (2.21) 

Just as in the L-T model, the set i? = in charged dust consists of the Big Bang/ Crunch 
singularity (which is now avoidable), and of the center of symmetry, which may or may 
not be singular. The conditions for the absence of a permanent central singularity are the 
following (not all of them are independent, but this is the full list): 

Nin) = = Q{r,)=Mir,) = M{n), (2.22) 
lim R/M^/'^ = (3{t) ^ 0, (2.23) 



lim r2(r) = 1 ^ lim E{r) = 0, (2.24) 

r^rc r^rc 

lim 2E/A^2/3 _ ijj^ (y2^^^ _ ^^2/3 drf ^ _ pQ^g^^ (2.25) 



and this last constant may be zero. 

In the mass-curvature coordinates, with A = 0, it is easy to solve the evolution equation 
dR/ds = u^. We quote here the only one case, for which we found in Ref. [1] that there 
may exist solutions avoiding both kinds of singularity, E < (see Fig. 5 in Paper I). 

We take i = +1 for expansion and i = —1 for collapse at the initial instant of evolution, 
and we denote 

2E ^ ' 

When < 0, a solution exists only with $ < M^/ (4£'^), and it is given by the parametric 
equations 



M /M2 ^ 
R = \ — ^ — cosu, 




6 



where is a parameter and is an arbitrary function of integration. This will avoid 

a BB/BC singularity only if M > and $ > 0. At s = sb, R starts off with the minimal 
value -Rmin = R+ (see eq. f l2.13p ). and periodically returns to this value, never going down 
to zero if -Rmin 7^ 0. The period is 

Tp = 27rM/(-2£;)3/2. (2.28) 

The period given by (12.281) is with respect to the proper time of the given shell of 
constant A4. To calculate the period in the time coordinate t, we must first transform 
the variables in eq. (12. 7p . to make {Ai,R) the independent variables and t the unknown 
function. We obtain 

rR- ^-CiM,R)/2^j^ 

- 1 + 2M/R + g2 (Q^^2 _ ^/^4) //^2 _ (i/3)Ai?2 



Note that each of the periods, Tp and T^, is a function of Ai only. We found in Paper I that 
even for a single bounce to be free of shell crossings, it is necessary that the bounce occurs 
in a time-symmetric manner, i.e. all the mass shells have to go through their minimal sizes 
at the same instant of the coordinate time t. 

Finally, we recall the 9 necessary conditions that the functions defining the model must 
obey in order that even a single nonsingular bounce can occur. We denote 

Fi 1 - (cVG) {Q,n' + QQ,NN ) , (2.30) 

and the conditions are: 

(1) E<0; 

(2) E > -1/2; 

(3) lim,_,^ F,/M'/' = 0; 

(4) r,^ < 0; 

(5) Q,N^ < Gjc^ at > and Q,jv^ = Gjc^ at = 0; 

(6) M = M- QQ,N^> 0; 

(7) M2 - 2EQ2 (g,^2 _ ^/^4) > 0; 

(8) 

W 

2E ^ ' 



(9) 

2E 



M + Jm^ - 2EQ^ _ ^/^4) 



^2.32) 



Conditions (1) - (9) must hold in a neighbourhood of the center. At the center, the left- 
hand sides of conditions (1) and (6 - 9) must have zero limits. In addition to this, all the 
regularity conditions at the center must be obeyed. Note that the conditions (I2.22p and 
(EMD, together with eq. (^TIB . imply that 



lim M/N = (2.33) 



this will prove useful later. 
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For practical calculations, the most convenient radial coordinate is N{r), and it will be 
used in some of the following sections. 

In Paper I we gave an example of a configuration that was supposed to obey conditions 
1-9 and all the regularity conditions. In fact, it did not obey the condition fl2.23l) . This 
oversight has already been corrected in an erratum, and in the version stored in the gr-qc 
archive, see Ref. [1]. That paper contained one more error: in the caption to Fig. 10 we 
stated that the two curves that seemed tangent were merely adjacent to each other. In 
truth, they were tangent - the minimal radius at bounce must be equal to the radius of 
the inner R-N horizon at every local extremum of the latter. The proof is given here in 
Appendix |Xl This is just a correction of an error that has no relation to the other results 
of the present paper. 



3 Inevitability of negative values of u 

From now on, we assume A = 0. 

From eq. fl2.1Up we see that if the quantity u defined in f l2.16p changes sign at a certain 
point (t, r), while R,r does not, then the energy density also changes sign at that point. We 
will discuss the consequences of the change of sign of u in the next section. In the present 
section we will show that, with the 9 conditions listed in Sec. 2 fulfilled, u is necessarily 
negative within an interval containing the instant of maximal compression (minimal size) 
of every mass shell on which Ad 0. We first prove that negative u appears in the 
hypersurface of maximal compression, and later we will show that u/R,r < for a certain 
period around the instant of maximal compression at every A4 ^ 0. 

At first sight, it seems that the conditions (1) - (9) given at the end of Sec. 2 secure 
M > in a vicinity of the center, since F — > 1, Q/N ±^/G/c^ and R/N^/^ /3 > as 
r — >■ Tc- However, the situation changes if i? = i.e. if one approaches the center along 
the locus of the inner turning points. Fig. 12 in Ref. [1] (and also Fig. H] further here) 
shows that along this locus R tends to zero faster than along any M. = const line, and eq. 
fl2.13p confirms this. At i? = i?+ we have 



u = r + 



(3.1) 



Q {Q,m' - G/c^) 

To guarantee e > 0, this should be positive everywhere, including the center = 0. This 
requirement can be easily fulfilled if Q,n'^ < G/c^ everywhere including the center. Then, 
it is enough to choose such Q{N) that the coefficient in front of the square brackets in 
f l3.1l) is small and tends to zero as A^ ^ 0, for example Q(A^) = ±(v^A'o/c2)e-"^', where 

a > is a constant and x '= N/Nq. 

In our case, when (5,iv niust obey condition (5), it turns out that the hmit of u at 
A^ — * is necessarily negative. This is seen as follows. We first observe that 

hm [-2EQ^ {Q,N^ - G/c")] /M^ = 0, (3.2) 

i.e. that the second term under the square root in (13. ip can be neglected compared to 
in the limit A^ — (see Appendix [B]) . Then, from (13.10 

2MQ,N 

lim u = 1 + lim — -. ^ — ^. 3.3 

^-0 ^-*^Q{Q,N^ -G/c^) ^ ^ 
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Since, by condition (5), Q^n (0) 7^ 0, we find 



lim u 





QQ,NN r — QQ,N r,Ar 



1 Q,n{0)T,n{0) 

2 2Q,jviv (0) 



Q,NN {o)-Q,N (o)r,Ar(o) 

2Q,NN (0) 



(3.4) 



(In deriving this, we first applied the de I'Hopital rule, then did an algebraic simplification 
in the result, and again applied the de I'Hopital rule in the denominator.) Now, by con- 
dition (5), if Q,N (0) > 0, then Q,nn (0) < 0, and if Q,n (0) < 0, then Q,nn (0) > 0. By 
condition (4), r,7v (0) < 0. Thus, if T,]\r (0) < 7^ Q,nn (0), then the second term in the 
last line of (13.41) can never be positive. If T,iy (0) < = Q,nn (0), then the limit in (13. 4p 
is —00 because Q,nn < in a vicinity of = 0. If r,iv (0) = 7^ Q,nn (0), then the limit 
is -1/2 < 0. If r,jv (0) = = Q,NN (0), then the limit of Q,n T,n /Q,nn at ^ is still 
negative because of conditions (4) and (5). Consequently, limTv^o"" < in every case, □. 

Since the inequality is sharp (actually, limAr^o"" < ~l/2), u will be negative already 
at some > 0. 

If we approach the center of symmetry along the locus of the outer turning points, 
R = then u remains positive up to the very center (see the proof in the final part 
of Appendix [B]) . This shows that the negative-w region does not exist permanently, but 
appears only for finite time intervals around the bounce instant. 

Thus, a nonsingular bounce of spherically symmetric charged dust is possible only when 
there exists a region of negative u for a finite time- interval around the bounce instant. We 
discuss the implications of this fact in the next section. 

4 Consequences of < 

The only place in the metric where u explicitly appears is via (12. 5p . which shows that 
e"^ = [R^r/uf- This is insensitive to the sign of u. One might thus suspect that it 
is enough to take e"^/^ = —R,r/u instead of (12. 5p to cure the problem. However, with 
such changed e^^^, the Einstein equations imply that C,r= —2(5(5, at e^/^/i?^ instead of 
(12. 6p (see the derivations in Refs. [1] and [1]), i.e. C is sensitive to the sign of u. The 
changes propagate through all the equations, and the resulting formula for energy-density, 
eq. (I2.10p . remains unchanged. This means that e becomes negative when m < 0, unless 
this change of sign of u can be offset by a change of sign of R,r or N,r. We show below 
that such an offset is impossible. 

The radial coordinate r is arbitrary, and all the formulae are covariant under the trans- 
formation r = f{r'), where / is an arbitrary function. The values of r label flow-lines of 
the charged dust, all flow lines of the same r form a 3-cylinder in spacetime, whose sections 
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of constant t are 2-spheres. Thus, for logical clarity, it is convenient to assume that the 
labeling is such that increasing r corresponds to receding from the center of symmetry^ 
Suppose that r = corresponds to the center of symmetry. Then, R,^ < is impossible 
in a vicinity of r = r^: since i? = at r = r^, R,r < at r > would imply i? < in a 
vicinity of the center, which is a geometrical nonsense. Thus, R,r > at the center, and 
so R,r /u < 0. But then the energy density will be negative at the center, unless N,r < 0. 
However, the physical interpretation of eq. (12. lip suggests that N,r cannot be negative: 
since = at the center, A^,,, < would imply < around the center. The final 
conclusion is that the energy density e/c^ must be negative where u < 0, at least in some 
neighbourhood of the center (but see further in this section). 

Is the set M = a singularity? The answer to this question depends not only on u, but 
also on the behaviour of R,r. If R,r = at m 7^ 0, then this is a shell crossing, which is 
a curvature singularity. The 9 conditions listed at the end of section [2] were derived from 
the requirement that R,r 7^ (actually, R,r /u ^ 0) throughout the evolution. All of them 
are necessary conditions for R,^ /u 7^ 0. A sufficient condition is not known, but anyway 
we wish to avoid R,r /u = 0, and will not consider this case here. 

If u/R,ry^ at M = 0, then there is no singularity at this location, but the energy 
density will be negative where u/R,r< 0. We will now investigate the behaviour of e in 
the set S'min, which is a 3-space composed of those points where R{t, r) assumes its minimal 
value i?+ given by (I2.13p . At those values, R,t = 0. 

We have already shown that u/R,r< (implying e < 0) in a vicinity of the center 
i? = in S'min- Where does S'min intersect m = 0? At the intersection, from (12.130 and 
(I2.16p . R must obey two equations: 



R = R+ = — ^ 
+ 2E 



N 



(4.1) 



r 

We substitute for M from (El]), recall that = 2E + 1, and solve (jH]) for Q,n- 

Q,N = ^[m-^/W^gqW) 

^R = Ri = M- - GQ^/c^. (4.2) 

This Ri is equal to the radius of the inner Reissner - Nordstrom event horizon corresponding 
to the mass shell 

Equation (14.21) should be understood as follows. The functions M.{r) and Q{r) are 
arbitrary functions in the model, and F is determined by M. via (12. 8p . Thus, (14. 2 p is an 
additional condition imposed on functions that are otherwise arbitrary, and the first of 
( 14. 2p may possibly have no solution. However, if ii has a solution, then this point must 
have the R value given by the second of (14. 2p . 



''This order of labels is reversed on the other side of a shell crossing, but we are considering models 
without shell crossings. 

^Formally, the solution of ()4.1|) includes also the outer R-N horizon, with "+" in front of the square 
root. However, here we consider eq. (|4.1|) within the set S'min that consists of the inner turning points of 
the various mass shells. Within this set, the equation R — M + \/ M"^ — GQ"^ /d^ has no solutions, as is 
known from the general properties of the R-N solution PJ. 
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Then we calculate i?,^ at the point determined by (14.21) . and obtain, using (12. 8p and 
the first of K2^ : 

iR,r = 0. (4.3) 

Thus, within the set S^i^ the locus of m = (if it exists at all) coincides with the locus 
of R,r = 0. However, the equations do not allow us to determine the sign of u/R,r at this 
location in the general case. In the explicit example given in Sec. 8, u/R,r < throughout 
S'min, but this may be a property of that particular model only. In general, we can only 
say that there exists a finite neighbourhood of the center R = within S'min in which 

U/R,r < 0. 

Analogously, we could consider the set S^ax consisting of those points where R = 
i.e. where all the mass shells attain their maximal radii. We showed in Appendix [B] that 
M > in a vicinity of the center within S'max- At the center, R,r > 0, so u/R,r > in a 
neighbourhood of the center. The whole calculation above can now be repeated with i?_ 
substituted for _R+, and the result will be that also within S'max the set m = coincides 
with R,r = 0, but this time at the value R = Ro = M. + a/A1^ — GQ"^ /d^. This means that 
there exists a finite neighbourhood of the center i? = within S'max in which u/R,r > 0, 
implying e > in the same neighbourhood. 

All this implies that a volume of charged dust evolves from negative energy density in 
S'min to positive energy density in S'max- Consequently, the regions of e < and e > are 
not separated by a comoving hypersurface - the matter particles proceed from one to the 
other. This should indicate that strong electric fields induce a hitherto unknown physical 
process inside matter. 

This concludes the discussion of the case when u/R,r^ at u = 0. 

Since we have found that u/R,r < in a neighbourhood of i? = in the hypersuface 
R = -R_|_ and u/R,r > in a neighbourhood of i? = in the hypersurface R = there 
must be some intermediate set between them on which u/R,r = = e. This set must be a 
subset of the m = set. Note that in a singularity- free model R,r can vanish only where 
u = but does not have to vanish everywhere on the m = set. Across the set on which 
u = ^ R,r, the energy density changes sign. In Sec. [8] we will trace all these boundaries 
on a numerical example. 

5 There is no geometric singularity at u/R,r= 

Now we consider the case when u/R,r= at u = 0. This is a singularity of the metric 
because \gu\ = |e^| — ^ oo there. However, none of the curvature or flow scalars become 
infinite at those locations, as we show below. 

The components of tensors referred to below are scalar components with respect to the 
orthonormal tetrad defined by the metric (12.11) : 

eO = e^/Mt, e^=e^/Mr, = Rd-d, = RsinMip. (5.1) 

These components are thus scalars. In this frame, the nonzero components of the Ricci 
^Because R,r = ^ u would be a shell crossing that was excluded by Conditions (1) - (9). 
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tensor and of the scalar curvature ar^ 



Rqo 
Rii 
R 



A + 



+ 



c4 R^R 



Rio — R. 



G 



4A ^ 



33 



2G uN,, 



uN, 



R'^ R^R. 



R^R' 



(5.2) 



The nonzero tetrad components of the Weyl tensor are 



Coioi 
C0202 

W 



dcf 



-C2323 = -2W^, 

Co303 = — C*1212 = 



Cl313 = 

_ ^ 



(5.3) 



None of these are singular at u/ R,r= 0. 

The scalars defined by the flow - the square of the acceleration vector u°'Ua, the ex- 
pansion 9 = u'^-^a and the shear a given by cx^ = ^a°'^aa/3, where 



are as follows 



N 



a 



2^3 



c/2 -A 
R 



(5.4) 

(5.5) 
(5.6) 

(5.7) 



The quantities ii^Ua and e-^/2^,i /i? are seen to be nonsingular at m = (the latter from 
(ETD), so we investigate e-^/M,t. We find from ([23]) and fl2Tl) : 



-C/2^ 



2e^/^ 

2e^ 
R,t 



2R,tr ^.QQ^N R, 
UR 



Q-^R.fR 



QQiN Rjr -C r> 2 
D2 ^ ^ 



it 



(^R 



R 



R 



u 

Rir 



M-QQ,Nr , Q^Q,n'-G/c^) 



i?3 



(5. 



^We substituted the solutions of the Einstein - MaxweU equations for C and A from (|2.5p - (12. 9p 
These formulae were calculated by the computer-algebra system Ortocartan [7] . 
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This is seen to be nonsingular at u/R,r= 0. Thus, all the flow scalars are nonsingular at 
u/ R^r = 0, just like all the curvature scalars. This shows that the singularity in gu = —e^ 
is merely a coordinate singularity. 
Note the following identity 

m'- A = (m^)' > 0. (5.9) 

Thus, at M = 0, A must be non-positive. This means, if the region m < exists, then the 
outer surface of this region either lies between the two R-N event horizons or touches one 
of them. (As shown before, the boundary of the region n < touches the horizons when 
u^ = atu = 0.) 

It is puzzling that this negative-w region is invisible in the mass-curvature coordinates. 
It is not surprising that the coordinate singularity atu = disappears - the mass-curvature 
coordinates turn thus out to be those that remove it. However, we showed above that at 
u/R,r = the mass density changes sign, and the change of sign of a scalar is an invariant 
property that should be visible in all coordinate systems. Equation fl2.20p shows that this 
can happen only by changing the sign of F. But then there is no way in which the e of 
f l2.20p could go through zero - F is finite everywhere, except at the turning points, at which 
Fu^ is finite. This seems to imply that the set m = is not covered by the mass-curvature 
coordinates. 

This is indeed likely. Equations fl2.5p - (12. 6p show that (12. 6p cannot be integrated 
across the u = set. Then, the definition of F via the comoving coordinates is F = 
R,r / (uu^Ai,r) , i.e. the transformation to the (A^, i?)-coordinates is singular at u/R,r = 
0. Thus, F should suffer a discontinuous jump from positive to negative values across 
the set u/R,r= 0, but this change remains invisible if we use the (A^, i?)-coordinates 
throughout. 

6 A vacuole around the center of symmetry is no rem- 
edy to u < 

Since, with the 9 conditions of Sec. [2], the region m < cannot be eliminated by a choice 
of the arbitrary functions, can we get rid of it by matching the charged dust metric to 
vacuum across a hypersurface r = const on the inside? i.e. by cutting an empty vacuole 
around the center of symmetry of the charged dust ball? We shall consider the matching 
only in the case A = 0. 

The electrovacuum spacetime matched to spherically symmetric charged dust must be 
Reissner - Nordstrom (R-N). However, if the R-N solution is to be applied around the 
center of symmetry, then it will contain a singularity unless its mass and charge parameters 
are both zero. In this case, the R-N solution reduces to the Minkowski spacetime. 

Can the charged dust metric be matched to the Minkowski metric, with all the matching 
conditions fulfilled? As eq. (12. lip implies, the matching requires that Q = A4 = at the 
matching sphere r = rj. Then, (12.70 shows that if r(rj) = 1, then R,t {t,ri) = 0, i.e. the 
charged dust particles that are initially on the matching hypersurface will remain on it 
permanently. Equation (12.60 shows that C,r{ri) = 0, thus C{t,ri) = Co(t) and can be 
made zero, so that e*" = 1, by a transformation of t. This all should happen at a nonzero 
R. 
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For the dust particles remaining permanently on the matching hypersurface, the mini- 
mal and maximal radius given by (12.131) should coincide and both be nonzero. In order to 
achieve this, the limit of —M/E at = must be positive, and in addition, to avoid -R+ 
bemg zero, the terms {M/2Ef and {Q,n^ - G/c^) /{2E) must be of the same order 
(so that their sum at r = rj does not equal {M/2EY). These conditions are impossible to 
fulfil - see Appendix O 

The conditions of a bounce at a nonzero R can be fulfilled in the cases > 0. However, 
then the arguments used in Refs. [21 H] and |1] still hold: even if Q = = occurs at 
i? > 0, the shell crossings will be unavoidable because the function F necessarily changes 
sign somewhere in the vicinity of the (5 = surface. 



7 A transient singularity at the center of symmetry 

Paper I and the present paper were concerned with avoiding the BB/BC, shell crossings 
and permanent central singularities. It turns out that there is one more kind of singularity 
that has slipped through all the tests applied so far. 

We noted in Sec. H] that i? — faster than anywhere else if we approach the center 
of symmetry along the hypersurface S'min; in which R = R^. The energy density given 
by (12.101) will be finite at the center provided the product uN,^ / {R^R,r ) < oo. We have 
already shown in Sec. [3] that u is finite at all points of the world line of the center, albeit 
negative when the center is approached from within S'min- Thus, to avoid a singularity 
at all points of the center, we must have lim^^re / {R^R,r ) = lim^^rc 3/ {R^) ,n < oo, 
i.e. lim^^rc (-R^) 7^ 0. However, as we show below, lim^^^^ (R^) ,Ar = if the center of 
symmetry is approached along the hypersurface S'min, and also lim^^^^ (R^) ,Ar = if the 
hypersurface S'min is approached along the center of symmetry. 

A nonzero limit of (R^) ,n at the center, where R = = N, means that R^ must tend 
to zero as fast as A^, i.e. that R must tend to zero as fast as N^^^ and no faster. Let us 
then take the first of (I2.27p . The sufficient condition for lim^^^^ R/N^/"^ 7^ is then 

M 

lim — ^ 0, (7.1) 

which does the job at those points where cos 00 7^ 1. With coscj = 1, however, where 
R = R+, we have 



R — R4. — 



2E) 



2EQ^ {Q,r,' - G/c^) 
M2 



(7.2) 



We proved in Appendix [B] that the fraction under the square root always tends to zero 
at r — > Tc if the other regularity conditions are obeyed. Thus, the expression in square 
brackets also tends to zero at r — > Tc, implying that R/N^^^ — > 0, i.e. a singularity at 

the center. 

We can approach the singular point along the world line of the center of symmetry. 
Assuming that lim^^^^ M/ {-2EN^/^) = /5o 7^ 0, we obtain 

Ym,^^^ = (3,{l-cosu)^Q (7.3) 
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(the dependence on time is now hidden in u). Taking the hmit of (17.31) as cosu 1, we 
again obtain the zero hmit of R/N^^^ at the feral point. 

Let us recaU: this transient singularity results because the i?+ given by (17. 2p goes to 
zero at r ^ Tc faster than N^^^. Can we avoid the singularity by forcing R^ to go to zero 
at r — s> as fast as N^^^ or slower? The following simple consideration provides a negative 
answer. Equation (17. 2p can be equivalently written as 



M 



2EQ^ {Q,n' - G/c') 
M2 



We can now use the reasoning in Appendix [B] to show that \imN-*o R+/N^^^ = indepen- 
dently of the forms of the functions Q, E and M, provided conditions (1) - (9) of sec. 2 
are obeyed. Equation flB.3l) tells us what forms p and (p,^^ — 1) must have in a vicinity of 
the center. Thus we have 

Q' {Q,n' - G/c') 
M 

x2 (1 + Apx'^"-' + (2ApX°--i + 



Ap (ap + l) x'^p + O^^ 
2r 2N 

-^ + 01^ +01. (7.5) 



Because of (17.51) and (17. ip . the expression under the square root in (17.40 has the limit 1 at 
r — > Tc. Consequently 

R: Ar2/3 

Note that since m < as the singularity is approached from within Smin, and at the 
center necessarily R > and R,n > 0, the density in the singularity becomes minus 
infinity. At the same time, if we start at a point of the center of symmetry different from 
the maximal compression instant, then, from (I2.2ip - (I2.24p and Q^N'^ijc) = G/c^ we 
have u{rc) = 1 > everywhere on A^ = 0, so e > everywhere on A^ = 0, including the 
limiting point of maximal compression. Consequently, e —>■ +oo at the point where the 
center of symmetry hits the maximal compression hypersurface. This singularity is thus 
direct ion- dep endent . 

This transient singularity exists within each hypersurface. However, with all other 
conditions of regularity obeyed, the limit of the period of oscillations of the solution (I2.27P 
at r — > Tc is infinite, as will be shown below. Thus, the transient singular point can be 
reached only once, where the bounce set s = sb{J^) hits the center of symmetry at a finite 
s. Every other hypersurface S^i^ escapes to the infinite future or the infinite past when it 
approaches the center of symmetry. 

Here is the proof that Tp — > oo, where Tp is the period of oscillation given by (12.280 . 

r^rc 

We know that M/{—2E) must tend to zero at the center as fast as A^^'''^. We know from 
(I2.33P that M tends to zero faster than A^, so let M oc A^-^+^, where e > 0. Consequently, 
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{—2E) tends to zero as fast as A^2/3+e_ This means that Tp = M/(— 2i?)^/^ behaves as 
N~^^'^, i.e. tends to infinity at the center. □ 

Since we were interested in sending a ball of dust through the Reissner - Nordstrom 
throat, we set up the initial conditions so that the state of maximal compression (minimal 
size) was attained simultaneously by all mass shells. In such a configuration, illustrated 
in Figs. 12 and 13 of Paper I (and also in Fig. H] further here), the hypersurfaces Smax 
(consisting of the instants of maximal size of all the shells) approach the center of symmetry 
in the infinite past and in the infinite future. Consequently, it is not in general guaranteed 
that the center of symmetry will ever reach the region where u > and e > 0, but in 
our explicit example, given in the next section, the region m < is contained in a small 
neighbourhood of S'min- 

However, we could take exactly the same example and set up the initial conditions so 
that the state of minimal compression (maximal size) is attained simultaneously. Then the 
dust ball would never emerge from the R-N throat because shell crossings would appear 
immediately after the bounce. However, initially and for some (perhaps infinite) time to 
the future and to the past, the center of symmetry would be in the region of positive u 
and e. 



8 An example 

To illustrate the statements of the preceding sections we will use the same example that 
we introduced in Paper I. To recall, in the example the function N was (and still will be) 
used as the radial coordinate. In the exemplary configuration, the charge function was 



Q{N) = p{x), (8.1) 

where q = ±1, to allow for any sign of the charge, x == N/Nq, Nq is an arbitrary constant, 
and 

Then 
and 



p{x) =x/{l + xy. (8.2) 



VG 1-x 



Fdx) = 1 - ^ {Q,n' + QQ,nn) = 1 - ^''[^^ly^ - (8.4) 

The original definition of the function E{x) was incorrect, as it implied a permanent 
central singularity. We quote here the corrected E{x) and the formulae dependent on it, 
from the gr-qc version of Paper I. The functions defined along the way obey conditions (1) 
- (9) of Sec. [2]- for proofs see the gr-qc version of Paper I. Thus 

bx^/^ 

where b is another arbitrary constant. With such E we have 

r(x) = , \ ^. , (8.6) 
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-0.2 



dp/dx 



2 4 6 8 10 

Figure 1: Graphs of the functions p{x), dp/dx, Fi{x), F2{x) and -^3(0;). The inset shows the functions 
p{x), F2{x) and ^{x) in a vicinity of x = 0. 



Mix) 



dx' 



def 



GNo 



fi{x). 



(8.7) 



Then, further: 



M 



Fo(x) 



G N 

M-QQ,Nr = —^F2ix) 



def 



x(l — x) 



(1 + x) Vl + bx^/^ 



For checking condition (7) we need the function F^lx) defined below. This condition is 
equivalent to 

F,{x) > 0, F,{x) = F2\x) - Fe{x), (8.9) 

where 



^^11/3 



(1 + te5/3) [l+xY 



1 - 



[l-x) 

:i+x) 



(8.10) 



Fig. [T] shows the graphs of the functions defined above with b = 2.5 (why this value - see 
below). This is a corrected version of Fig. 7 from Paper I. 

For checking the remaining conditions we need the functions -^4(3;) and F^lx) defined 
below. Condition (8) is equivalent to: 



^4(0;) > 0, where 
and condition (9) is equivalent to: 



FJx) = FAx) 



6xVl + bx^/^ 



Foix 



(8.11] 



F^lx) > 0, where F^lx) =^ Fi{x 



6xVl + bx^l'^ 
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F2{x) + ^/F^]. (8.12) 




0.5 1 1.5 2 



Figure 2: Graphs of the functions ^4(2;) and Fz{x) with 6 = 2.5. The graphs of Fi{x), F2{x) and ^3(2;) 
are shown for comparison and for scale. 



Fig. [2] shows the graphs of F4{x) and F^i^x) with b = 2.5. 

The reason for choosing 6 = 2.5 was this: the graph of F5 is sensitive to the value of b. 
With decreasing b, the local minimum of F5 becomes smaller, and for b small enough (for 
example, 6 = 0.75) F5 < around the minimum. With b > 2.5, the minimum is clearly 
positive. The value b = 2.5 was used in all the subsequent figures. 

Figs. 9 and 10 of Paper I do not change in any noticeable way after the correction done 
in E{x), so they need not be repeated here (see the gr-qc version of Paper I). However, 
Figs. 11 and 12 do change significantly (also in consequence of the different value of b 
adopted here), and their correct versions are shown in Figs. [3] and H] here. Fig. H] also 
shows the light cone of the transient singularity at {t,R) = (^^,0) - to demonstrate that 
all mass shells enter this cone soon after going through the minimal size. 

The other qualitative conclusions presented in Paper I remain unchanged. Below, we 
present the new results found in the present paper. In all the subsequent figures, the 
radial coordinate r was defined as r = N^^^, so that one can instantly see whether there 
is a central singularity {R,r = at the center) or not {R,r > at the center). Fig. [5] shows 
the functions u{x) and R,r (x) (left graph) and the energy density e(x) (right graph) in 
the surface R = R^. The functions u{x) and R,r (x) have a zero at the same value of 
X and opposite signs in all other points, except at x = where R,r and m < 0, so 
u/R,r — > —00. As follows from the calculations in Sees. |l]and[7l in this surface the energy 
density is negative in a vicinity of the center and goes to —00 at the center. 

The fact that the point where u = = R,r does not seem to be singular is rather 
mysterious because it is a limiting point of the contour R,r = in the (t, x)-plane, and the 
rest of the contour is a shell crossing (see Fig. [7]), where the energy density does go to 
±cx). But the curves in Fig. H] are all for smaller values of x, so the singularity is not in 
its range. 

Fig. [6] shows the functions u{x) and R,r (x) (left graph) and the energy density e(x) 
(right graph) in the surface R = -R_. Within this surface u{x) and R,r (x) are everywhere 
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Figure 4: The curves R{Ai,t) corresponding to several values of M. The mass increases uniformly from 
X = 0.01 on the lowest curve to a; = 0.1 on the highest curve. The origin t = R ~ is the point singularity 
discussed in Sec. [T] A null cone with the vertex at the point singularity is also shown - as can be seen, 
by the time the outermost mass shell reaches its maximal size, it is within the future light cone of the 
singularity. The inset shows the light cone in a small neighbourhood of its vertex - all mass shells enter 
this cone soon after going through the bounce. The region u < 0, as follows by comparison with Fig. [71 
would occupy a very small neighbourhood of the origin and would not be visible at the scale of this figure. 
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Figure 5: The functions u{x) and R,r (x) (left graph) and energy density e{x) (right graph) in the surface 
R = The first two change their signs at the same value of x, which causes that u/R,r< in this 
whole surface, implying e < in the whole surface. The inset in the left graph shows the neigbourhood of 
the intersection point. Note that R,r^ at r ^ 0, which causes the singularity at the center. As shown 
in Sees. U and [71 the energy density in this surface is negative in a vicinity of the center and goes to — cxd 
at the center. 

positivelfl The energy density is everywhere positive and finite (it does not tend to zero, 
but to a small positive value as r ^ 0, as shown in the inset). 

Fig. [7] shows the curves u = (lower part) and R,r = (upper part) in the (t,a;)-plane 
(actually, each curve is one half of the full contour, which is mirror-symmetric with respect 
to the X-axis). The function u is negative to the left of the lower contour and positive to 
the right of it. The function R,r is negative to the left of the upper contour and positive 
to the right of it. As stated above, the set R,r = is a shell crossing. Its presence proves 
that the 9 necessary conditions listed at the end of Sec. [2] were not sufficient. Note the 
characteristic features of the curves, consistent with the calculations of Sees. H] and O 

1. The lower curve hits the t-axis (the center of symmetry) at t = (in the figure, 
we chose = 0, so t = is the simultaneously achieved state of minimal size). Thus u 
remains positive all the time as we proceed toward t = along the center of symmetry. 

2. This curve intersects the x-axis at some a; > 0. This intersection point coincides 
with the point where the u{x) and R,r (x) curves intersect in Fig. [51 The upper curve 
begins at the same point. 

3. To the right of the lower curve, we have u > 0, and to the right of the upper curve 
R,r > 0. The energy density is positive in the area to the right of both curves, and negative 
to the left of any of them, including the x-axis. 

In order to better visualise the variation of the functions u, R,r and e in the {t, R)- 
plane, we provide three further graphs. Fig. [H] shows the functions u{x) and R,r (x) along 
the line t = 0.004 in Fig. [71 Fig. [9] shows the energy density along the same section. It 
changes sign every time when the line t = 0.004 crosses one contour or the other. Where 
M = 0, it changes from negative to positive by smoothly going through zero; where R,r = 
it changes sign by jumping from +oo to —oo or the other way round. 

^We recall that the surface R = R- in our example never intersects the worldline of the center of 
symmetry, but approaches it asymptotically a,s t +oo or i — > — oo, depending on whether it lies to the 
future or to the past of i? = R+. 
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Figure 6: The functions u{x) and R,r (x) (left graph) and energy density e{x) (right graph) in the 
surface R = R-. The functions u{x) and R,r (x) are positive everywhere, implying e > in the whole 
surface. This time R,r > at r = 0, i.e. there is no singularity at the center - as confirmed in the right 
graph. The inset in the right graph shows that the energy density at the center is not zero, but has a 
small positive value. 
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Figure 7: The u = and R,r = curves in the {t, x) plane. Each of the functions u and R,r is negative 

to the left of its graph and positive everywhere to the right. The full contour is mirror-symmetric with 
respect to the x-axis. The R,r = curve is the locus of shell crossings. The inset shows a closeup view 
of the neighbourhood of the origin - to demonstrate that the u = curve hits the f-axis only at t = 0. 
Thus, at the center of symmetry, u > everywhere except at the singular origin (where its value becomes 
dependent on the path of approach). See more explanation in the text. 



21 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.002 0.004 0.006 0.008 0.01 

X X 

Figure 8: The functions u{x) and R,r (x) along the hne t = 0.004 in Fig. \7\ The left graph presents 
these functions for larger values of x, the graph on the right for smaller values of x. Note that u increases 
to a positive value at the center x — 0, even though it is negative in a vicinity of the center. 




Figure 9: The energy density e{x) along the line t — 0.004 in Fig. [T] The left graph presents e{x) 
for large values of x, the graph on the right for small values of x. Note, in the left graph, that the curve 
goes through zero as it crosses the u — contour, but very near to that point it crosses the R,r = 
contour which is a shell crossing. There, it goes to +oo on one side and to —oo on the other. The second 
intersection with R,r = is also seen at x ~ 2, where the density changes sign again by a jump from 
— cxD to +0O. The right graph shows how the density becomes negative with very large absolute value in a 
vicinity of the center at a; = 0, but then, still closer to the center, grows to very large positive values. 
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Figure 10: The functions u(x) and R,r (x) (left graph) and the energy density e{x) (right graph) along 
the hne t = 0.025 in Fig. [7l 
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Figure 11: The functions u{x), R,r (x) and e{x) along the line t = 0.05 in Fig. [71 



Fig. [To] shows the same functions along the line t = t2, where ^2 hes between t = 0.004 
and the right end of the the u = contour in Fig. [3 t2 = 0.025. Finally, Fig. [TT] shows 
the functions u{x), R,r {x) and e(x) along the line t = t^, where lies completely to the 
right of the u = contour in Fig. [71 In each case, the functions behave exactly in the way 
in which Fig. [7] implies they should. 

Fig. 13 of Paper I, which showed the schematic Penrose diagram of the evolution of our 
exemplary model, is still qualitatively correct. However, now it has to be supplemented 
with a few new elements - the images of the surfaces u = and the light cone of the point 
singularity at the center. The complete version of that diagram is shown in Fig. [121 and 
here is the explanation. 

The diagram is written into the background of the Penrose diagram for the maximally 
extended Reissner-Nordstrom spacetime (thin lines). C is the center of symmetry, Sb is the 
surface of the charged ball, Srn is the Reissner-Nordstrom singularity. The interior of the 
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Figure 12: A schematic Penrose diagram for the configuration defined by eqs. ()8.1|) - (|8.2p and (|8.5p . 
The image of the interior of the dust ball is deformed to make its center R = coincide, wherever possible, 
with the set r = of the background Reissner - Nordstrom diagram - this is why the null geodesies there 
are not straight lines at 45°. See more explanation in the text. 
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body is encompassed by the lines C, E, Sb and B; the only singularity that occurs within 
this area is the point singularity at the center, marked PS. The dotted lines issuing from 
PS mark the approximate position of the future and past light cone of the singularity (i.e. 
the Cauchy horizon of the nonsingular part of the surface R = i?+ ) . Lines B and E connect 
the points in spacetime where the shell crossings occur at different mass shells. Nl (N2) 
are the past- (future-) directed null geodesies emanating from the points in which the shell 
crossings reach the surface of the body (compare Fig. H]). The line Sb should be identified 
with the uppermost curve in Fig. |H The top end of Sb is where the corresponding curve 
in Fig. m first crosses another curve, the middle point of Sb is at t = in Fig. |H The two 
dotted arcs marked "m" symbolise the hypersurfaces e = 0, on which the energy density 
changes from positive (below the lower curve) to negative and then to positive again (above 
the upper curve). The distance between these arcs is greatly exaggerated; at the proper 
scale of the figure they would seem to coincide. The arcs are parts of the u = contour 
shown in Fig. [71 In Figs. H] and [T^l the surface of the body is at smaller x than the top 
of the u = contour in Fig. [71 Every 7^ world line of the dust necessarily enters the 
region of negative energy density for a finite time interval around the instant of maximal 
compression. 



9 Permanent avoidance of singularity is impossible 

In Paper I we hypothesised that a permanently nonsingular configuration of charged dust 
might exist - provided the period of oscillations is independent of the mass Ai. We show 
below that in our class of configurations the period cannot be independent of M. because 
conditions (1) - (9) prohibit this. 

The condition Tc^m = 0, if considered by explicitly differentiating (12.291) . leads to a 
very complicated integral equation that we were not able to handle. (The integral in 
fl2.29p cannot be explicitly calculated without knowing the explicit form of C{Ai,R).) 
However, let us recall that both Tp given by fl2.28p and Tc given by fl2.29p are in general 
functions of Ai only. Thus, we conclude that Tc is a function of Tp, and so the condition 
for the special case of Tc being constant is 



Hence, will be constant when 



dM dTp dM 



dTj, dT„ 

-0^—7 = 0, (9.2) 



dM dN 

since f[2:2TD implies that A^,^ ^ 0. From I^M), <B) and Tp = C = const we then find 

(_2E)3/2 = (1 - r^) = C^/^ {M - QQ,N r) . (9.3) 

From here 

1 /i-r2\'/2 M 



QQ,^=--i^—^j +—■ (9.4) 
Given r(A^) and Ai{N) found from (12.211) . this determines Q{N) simply by integration. 
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We now verify whether (19.41) is consistent with the regularity conditions (12.221) - (I2.25P 
and the 9 conditions hsted after (12.301) . 

From (12.281) we see that Tp = C = const imphes 

Um ^ = (9.5) 

where, as stated after (12. 5p . D may be zero. Then, comparing (12.231) with (I2.27p . we see 
that 

M , , , 

ii^o (-2E)A<V3 ="^Q' "<^' ^9-6) 

where a is a constant. This, together with (12.251) . means 

M 

hm — — = a < oo, (9.7) 

i.e. that D ^ But with CD ^ 0, eq. (19.51) is in contradiction with \im.M_,QM/M ^ 
(see the remark at the end of Sec. [2]). 

We have thus proven that a configuration with E <{] cannot pulsate with the period of 
pulsations being independent of the mass M., while being singularity- free for ever. Different 
periods for different mass shells will necessarily cause shell crossings, during the second 
cycle at the latest (Fig. H] illustrates this). 

Permanently nonsingular oscillations of weakly charged dust, with the period indepen- 
dent of mass, are possible only when there is a central singularity. An example of such a 
configuration results when F = 1/ cosh(6x), where 6 is a constant. 

A spherically symmetric charged dust configuration can be permanently nonsingular 
only if it is static. Such configurations do indeed exist with special forms of the arbitrary 
functions, as pointed out in Refs [1] and [1], but they are not interesting from the point of 
view of avoiding singularities. 



10 Conclusions 

The conclusion of this paper is: the weakly charged spherically symmetric dust distribution 
considered here {Q,n'^ < G/c^ at > and Q,n'^ = G/c^ at N = 0) must contain at least 
one of the following features: 

1. A Big Bang/Big Crunch singularity; 

2. A permanent central singularity; 

3. A shell crossing singularity in a vicinity of the center 

4. A finite time interval around the bounce instant in which the energy density becomes 
negative, and a transient momentary singularity of infinite energy density at a single point 
on the world line of the center of symmetry. 

A fully nonsingular bounce is possible for a strongly charged configuration, Q,n'^ > 
G/c^, which can exist only with E > Q - see the comments after (I2.14p . This would be a 
collapse followed by a single bounce and re-expansion to infinite size. This phenomenon 
occurs also in the Newtonian limit - the bounce here is caused by the prevalence of elec- 
trostatic repulsion over gravitational attraction. Such an example was reportedly found, 
but never published, by Ori and coworkers [6J. 
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A permanently nonsingular pulsating configuration of spherically symmetric charged 
dust does not exist. At most, a single full cycle of nonsingular bounce can occur, and 
shell crossings will necessarily appear during the second collapse phase. The nonsingular 
bounce occurs at i? > 0, but the momentary isolated singularity at the center of symmetry 
is still there. 

The possibility of e going negative in the presence of electric charges does not seem to 
have been noticed and may need further work on its interpretation. If the negative energy 
density region existed permanently in some part of the space, with comoving boundary, 
then we might suspect that this is a consequence of a bad choice of parameters that implies 
unphysical properties in that region. However, here we have the case in which the energy 
density is positive for some time, and then these same matter particles acquire negative 
energy density in a time-interval around the bounce instant. This suggests that there is 
some physical process involved in this, which should be further investigated. 

A question arises now. The uncharged limit of the family of configurations defined by 
eqs. diH) - flCT]) . Q = 0, is the Lemaitre - Tolman (L-T) model PE]. For the latter, 
shell crossings can be avoided, as is well-known since long ago [10]. Why, then, are they 
unavoidable with Q 7^ 0? 

The answer is this: in the L-T model, the Big Bang (BB) or Big Crunch (BC) are 
unavoidable (both are present when E < 0). In the cases that are colloquially called 
"free of shell crossings" , in reality the shell crossings are not removed, but shifted to the 
epoch before the BB or after the BC, or both. Thus, the shell crossings are no longer 
in the domain of physical applicability of the model. When the BB/BC singularities are 
replaced with a smooth bounce in charged dust, the shell crossings that were hidden on 
the other side of BB/BC become physically accessible, and they terminate the evolution 
of the configuration. 



A The proof that R^{N) coincides with H^i^{N) at the 
local extrema of H^i-^{N) 

We stated in the caption to Fig. 10 in Paper I, that the curves ifmin(A^) and R^{N) 
were not really tangent at the point where ifmin(A^) has its maximum, but just close to 
one another. In fact, they not only were tangent at that point, but had to be tangent, 
independently of the explicit forms of the two functions. We show here that this is a 
general law: at every local extremum of H^i^{N) (call it iVg) we have ifinin(A'e) = R+{Ne). 
The radius of the inner R-N horizon as a function of N is 

^^min(A^) =M- ^M^ - GQyc^ (A.l) 

where Ai and Q are assumed to depend on through x = N/Nq. We do not assume any 
explicit form of the functions and Q, we only use the general properties (12.91) . (I2.15p 
and fl2:23|l . 

At every local extremum N = N^. we have H^i^ f^f = 0. Using (12.231) . this means 

rVA12 _ GQ^cA =[MT-QQ,nI, (A.2) 
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where the subscript e denotes the value at the extremum. From here 

{2MT - QQ,N 



'GQ^' 




QQ,N 




e 


Y2 



Substituting this value of GQ^/c^ in flA.ip wet get 

H^UNe) = QQ,N /r. 

Substituting the same value of GQ^/c^ in (12.131) for i?+ we find 

R4N,)=QQ,N/T = H,^UNe). 

□ 



(A.3) 



(A.4) 



(A.5) 



B Proof of ( 13:21 ) 



For simplicity, to avoid physical coefficients, we introduce the functions p{x), and 
F2{x) (where x = N/Nq, Nq = const) by 



^ / „ r\ dcf 

Q{N) = ± — p{x) 



def GNq 
M = fi{x) 



M{N) ^ F,{x) 



Then, to prove (13.21) we must prove that 



lim [-2V [pJ - 1)] /i^2' = 0. 
By condition (5), in a vicinity of a; = 0, p must be of the form 



p{x) 



X + Apx"'' + a 



- 1 = 2ApX°--^ + C)„^_i 



(B.l) 



(B.2) 



:b.3) 



where Ap and > 1 are constants, and 0^^ is an unspecified function with the property 
lima;_>o Oq-p /x"*" = 0. Then, by the regularity condition (12.2 



r = 1 + Aex""^ + O, 



-2E = -2Aex''^ + 0, 



0B5 



where > 2/3, by ([225D. It follows from fICT]) and (^M> - O that 

fi{x) = X + A^x'^^^' + Oa^. 

Hence, 

F2 = fi- pp„ T = {A^- Ae) - Ap («p + 1) x"- - - 



:b.4) 



(B.5) 



:b.6) 



Thus F2 is of order ap or aE + 1 > 5/3, whichever limit is lower. However, the limit via 
aE leads to a central singularity, since then, from (12.271) . R/N^^^ — >• as ^ 0, so (I2.23P 
will not hold. Consequently, F2 is of order ap, while the second term under the square root 
in (13.11) . —2Ep^ {p,x'^ — 1) is of order 0^2 = aE + ctp + 1- Suppose that 0^2 < 2ap, so that 
the second term under the square root in (13. ip is of lower order than the first one. This 
translates to > a^; + 1. But condition (I2.23p . in combination with (12.270 . requires that 
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oip = aE + 1/3, or else a central singularity will appear. Thus, —2Ep'^ {p,^^ — 1) must be 
of higher order in x than in eq. (13. ip . i.e. (13. 2p must hold. □ 

Now consider the case when we approach the center of symmetry along the locus of 
the outer turning points, R = Then m_ = F — QQ,n I R- remains positive up to the 
center, as is seen from the reasoning below. The conclusion that —2Ep'^ {Pix^ ~ 1) is of 
higher order in x than still applies, so we have 

2EQQ 

= 1 + Q,^(0) hmEQ/M= l + Q,;v(0)limEp/F2. (B.7) 

We know from the above that p{x) = x+ApX'^p+Oap, F2 = —Ap (op + 1) x'^^ — Oa^ and E = 
Aex°''^ + OaE, while the regularity condition (I2.23P implies, via (12.270 . that ue = a^ — 1/3. 
This, taken together, implies that limx^o Ep/ F2 = 0, i.e. that limAr^o^- = 1 > 0, i.e. 
that U- remains positive up to the very center. 



lim U- = lim 



C Impossibility of a vacuole around R = Q with ^ < 

Let us assume E{x) as in (1B.4P and p{x) in a form similar to (1B.3P : 



p{x) = ax + Apx"p + Oa^,^ p,^ = a + Apapx"^ ^ + Ca^-i, (C.l) 
For the function fi{x) this implies 

fJ^ = x + -^^x''^^' + (C.2) 
aE + I 

Then 

F2{x) = ^i-pp,^T = il-a^)x- Ae [20^ ^ ) - aAp (a^ + 1) x"'' 

\ «£; + 1 / 

+ O^^ + O^^+i. (C.3) 

Now we have two cases: o? ^ \ and o? = \. In the first case we have F2 = (1 — a^) x + 
0\{x). Thus, in order that MjE has a nonzero limit at x — > 0, it follows that a^; = 1. 
Then, if the terms {M/2Ef and (q^^2 _ ^^^4j i^2E) in fl^Tgj) are to be of the same 
order, the ratio p^ {p,x'^ ~ 1) / should have a nonzero finite limit at x ^ 0. But, with 
o? ^\ and as = 1 we have 

{pj - 1) _ x^ (g + Apx^p-^ + (g^ - 1 + 2ApapX'^p-^ + 

2E 2Aex + d 

0, (C.4) 



which has the consequence that -R+ = 0, i.e. the conditions Q = = cannot be imposed 
at > 0. 

In the second case, = 1, we can assume a = 1 (nothing in the equations depends 
on the sign of Q); then F2 is of the order of either x"^"*"^ or x"'', whichever exponent is 
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smaller. With (a^; + 1) being smaller, the result is immediately seen: If F2 is of the order 
^OB+i ig of the order then limx^o M/E = 0, and lim2;^o-R+ = 0. With ap 

being smaller we have 

p2 (pj - 1) x^^+i (1 + Apx'^"-' + O^^.iY {2Apap + Oo) 



2E 2Aex<^e + O, 



(C.5) 



CUE 



This will have a finite limit at a; = if a^; = + 1. But we are considering the case when 
F2 is of order x"'', which means that the limit of F2/E, i.e. the limit of M/E at a; — > will 
be infinite, which again shows that the matching conditions cannot be fulfilled at a finite 
nonzero R. □ 

Acknowledgement We are grateful to the referee for pointing out several errors in an 
earlier version of this paper. 
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